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We study bosonic atoms with two internal states in artificial gauge potentials whose strengths are 
different for the two components. A series of topological phases for such systems is proposed using 
the composite fermion theory and the parton construction. It is found in exact diagonalization that 
some of the proposed states may be realized for simple contact interaction between bosons. The 
ground states and low-energy excitations of these states are modeled using trial wave functions. The 
effective field theories for these states are also constructed and reveal some interesting properties. 


I. INTRODUCTION 

Atomic gases in the quantum degenerate regime pro¬ 
vide unprecedented opportunities for studying quantum 
many-body systema^S. These charge neutral atoms do 
not couple to electromagnetic fields as charged particles 
do. To this end, various methods for generating arti¬ 
ficial gauge potentials have been proposed and some of 
them have been successfully implementecP*^. A concep¬ 
tually simple way of generating an effective magnetic field 
for neutral atoms is using the Coriolis force in a rotat¬ 
ing systerrP^. There has also been substantial progress 
in generating Abelian and non-Abelian artificial gauge 
potentials usirm atom-light coupling in continuum or on 
optical lattice^®!^. 

The synthesis of artificial gauge potentials extend the 
realm of phenomena that can be simulated using cold 
atoms to topological phases of matte r. One primary ex¬ 
ample is the quantum Hall stateJ^^E^ which occur in two- 
dimensional electron gases exposed to an external mag¬ 
netic field (i.e. an Abelian gauge potential). Another ex¬ 
ample of interest in rece nt yea rs is time reversal symmet¬ 
ric topological insulator j^^b^ l arising from spin-orbit cou¬ 
plings (which can be interpreted as a non-Abelian gauge 
potential). The cold atom systems allow us to study not 
only fermions but also bosons. The interaction potential 
between bosons can be modeled very well as the short- 
range contact interaction in most cases, but long-range 
interactions can be found in dipolar moleculeJ^. 

The topological phases of bosons in artificial gauge 
potentials have been studied in many previous works. 
For one-component bosons, composite fermion states 
and Moore-Read state can be realized using contact 
interactiorPsHsD while Read-Rezayi stat es m ay appear 
if long-range interactions are introducecP^EH Pqj. |;.y^ro- 
component bosons, composite fermion states can also 
be realized using contact interactiorP^I. At filling fac¬ 
tor V = 4/3, a non-Abelian s pin-sin glet state has been 
studied as a possible candidat^^SUl but there is a com¬ 
peting composite fermion stat^^. In addition to frac¬ 
tional quantum Hall (FQH) states, a symmetry pro¬ 
tected topological (SPT) stat e with out anyonic excita¬ 
tions can be realized at jz = j 24 | 27| l29] interesting 


possibility for two-component bosons is creating artifi¬ 
cial gauge potentials that depend on the atomic species; 
e.g., the direction of th e pot ential for an atom is related 
to its magnetic momenlPIBl. This is partially driven by 
the physics of two-dimensional topological insulators, for 
which a minimal model consisting of integer quantum 
Hall (IQH) states of spin-up electrons in positive mag¬ 
netic field and spin-down electrons in negative magnetic 
field can be constructed. The quantum phases of two- 
component bosons in equal b ut op posite synthetic mag¬ 
netic fields have been studiecP^lMJ jn general, one may 
ask what could happen when multi-component bosons 
are coupled to general non-Abelian gauge potentials. In 
view of the successful identification of various topolog¬ 
ical phases in one- and two-component bosons, we can 
expect to see a plethora of exotic phenomena including 
FQH and SPT states due to the interplay between in¬ 
ternal degrees of freedom, artificial gauge potentials, and 
tunable interactions. Theoretical investigations of such 
systems would be useful guidance for future experimental 
studies as they can help to narrow down the parameter 
space one should explore. 


In this paper, we consider the cases where two types 
of atoms experience artificial gauge potentials pointing 
to the same direction but having different magnitudes. 
A straightforward way of creating such a system is to 
use a superposition of two potentials with different mag¬ 
nitudes: one has the same direction for both compo¬ 
nents and the other has opposite directions for the two 
components. T he co mposite fermion theorjP^ and par- 
ton constructioiP^M! will be used to construct trial wave 
functions describing topological phases in such systems. 
These trial wave functions are compared with exact di¬ 
agonalization results. The physical properties of these 
states are also studied using effective field theories. In 
some cases, the systems can be interpreted as “topolog¬ 
ical insulators” of emergent particles in the sense that 
the effective magnetic fields experienced by two types of 
composite fermions have the same magnitude but point 
to opposite directions. 
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II. MODELS AND METHODS 

Although some experimental systems implement arti¬ 
ficial gauge potentials in optical lattices, we shall con¬ 
sider continuum models for simplicity. This is a reason¬ 
able and convenient approximation because the Harper- 
Hofstadter model for particles moving in a periodic po¬ 
tential with magnetic fluxe^^lHlIl reduces to a continuum 
model when the flux per plaquette is small enougfP^. 
For a system of bosons with two possible internal states 
(called spin-up and spin-down) in two dimensions, the 
single-particle Hamiltonian is 


only non-vanishing component is the zeroth one Vq. We 
choose ga-T = 9 = 47r£| such that Vq = 1 in the spin-up 
wave function basis and it is used as the unit of energy in 
what follows. It is assumed that the bosons are confined 
to their respective LLLs and mixing with higher Landau 
levels is negligible. In the second quantized notation, the 
many-body Hamiltonian is 

E (4) 

<7T {mi} 

where the coefficients can be calculated using 

the wave functions (see Appendix A for details). 


_ I (p —e^A'f)^ 0 

2M 0 (p —e^Af)^ 


( 1 ) 


A. Sphere 


where e'^A®’ is the artificial gauge potential for bosons 
with spin a and M is their mass. The introduction of 
“charge” e'^ for the particles is redundant at this stage 
but will be useful later when we discuss the physical inter¬ 
pretation of our results. The physical quantities should 
only depend on the product e'^A°’. The magnetic field 
strength is V x A®’ = We will focus on 

systems in which and have the same sign 

but different magnitudes {e^B^ is always taken to be 
smaller). The solution of the single-particle problem is 
two sets of Landau levels (LLs) for spin-up and spin- 
down bosons, respectively. In the symmetric gauge with 
A°' = H‘^(—2//2,a:/2,0), the lowest Landau level (LLL) 
wave functions are 




z™exp 

C’^V27r2'"TO! 


( 2 ) 


where z = x + iy is the two-dimensional complex coordi¬ 
nate, m is the z component of the angular momentum, 
and £„ = y/li/(e‘^H‘^) is the magnetic length. In later 
discussions, we will consider systems in which the parti¬ 
cles carry negative charge or see a magnetic field along 
the —Cz direction. Their wave functions can be obtained 
by taking the complex conjugate of Eq. 

For a many-body system with bosons in the spin 
state (T, we perform exact diagonalization to find its 
ground state and low-lying excitations. The bosons are 
placed on a compact manifold such as sphere or torud^QUil 
to eliminate edge effects. The number of magnetic fluxes 
through the surface of sphere or torus for the spin state 
a is denoted as N^. The filling factors for the two 

spin states are defined separately because The 

magnetic lengths for the two spin states are differ¬ 
ent and we choose the smaller one £{- as the length scale. 
The two-body contact interaction potential between the 
bosons is 


V(ri-ri} = Sir.-ri)[fyn) (3) 

where g^-r’s characterize the interaction strengths. In the 
pseudopotential representation of this interactiorl^, the 


For the spherical geometry, a radial magnetic field can 
be generated by a magnetic monopole at the center of 
the sphere. The LLL wave functions on the sphere are 
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+ Y ] 
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where u = cos(0/2)e*^/^, u = sin(0/2)e“*^/^ are spinor 
coordinates {9 and ^ are the azimuthal and radial an¬ 
gles in the spherical coordinate system) and m is the z 
component of the angular momentum. The radius R of 

the sphere is £^ = £^ A"^/2. The matrix element 

Fmlm 2 m 4 m 3 Smi+m 2 ,m 3 +m 4 where 6ij is the usual Kro- 
necker delta function. This means that the z component 
of the total angular momentum is a conserved quantity. 
The many-body eigenstates can be further classified by 
their total angular momentum eigenvalue L{L + 1). 


B. Torus 


For the torus geometry, we consider a rectangle torus 
spanned by the basis vectors Li = Liex and L 2 = L 2 ey 
and choose the vector potential A*^ = (0, 5*^2;, 0). The 
LLL wave functions on the torus are 


ijj 


NZ 


m 




(to -I- kN^) 


1 

2 


X 

£b 



( 6 ) 


where the magnetic length £„ = ^L 1 L 2 /The 

matrix element cx where S^ J 

is a generalized Kronecker delta function defined as 

S^j = 1 iff z mod = j mod (7) 

with being the greatest common divisor of and 

This means that the many-body eigenstates can be 
labeled by a special momentum quantum number K = 

SE K) mod Nf. 
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III. TRIAL WAVE FUNCTIONS 


As the parameter space of our problem is quite large, 
we turn to the composite fermion theorjEH for guidance 
about where to search for gapped topological phases. 
This theory was originally developed for electronic sys¬ 
tems, where one electron binds an even number of mag¬ 
netic flux quanta (one flux quantum is (j)o = hc/e) to 
become one composite fermion. The composite fermions 
move in effective magnetic field and form IQH states in 
many situations because they generally interact weakly 
with each other. For bosons in artificial gauge potentials, 
it is not immediately clear how to define “flux quantum” 
as they do not carry charges. The simple choice used 
in previous works is to take the charge of the bosons 
to be 1 and attach an odd number of magnetic flux 
quanta to each boson. For one-component bosons and 
two-component bosons in the same magnetic fluxes, this 
prescription yields accurate appr oxima tions to the exact 
eigenstates of contact interactiorPlI^. If this is still ap¬ 
plied when the two components experience different ar¬ 
tificial gauge potentials, the magnetic field strengths for 
spin-up and spin-down bosons would be different. For 
certain system parameters, we will have cases in which 
the effective magnetic helds for spin-up and spin-down 
composite fermions point to opposite directions. This 
can not happen if the two components experience the 
same magnetic fluxes as studied befor^^. 

The heuristic picture used above is captured by the 
trial wave functions 




t 

^/t 



Xi 


( 8 ) 


where Xf is the wave function of / completely filled LLs 
(without the Gaussian factor). The X/.^ dX/j^]*) fac¬ 
tor only contains the coordinates of spin-up (spin-down) 
bosons, which describes spin-up (spin-down) composite 
fermions in positive (negative) effective magnetic field. 
The factor xi = Y[j<ki^j ~ contains the coordi¬ 
nates of bosons in both spin states and implements vor¬ 
tex attachment for all of them. The reason for using a 
sim sign rather than an equal sign is two fold. First of 
all, we do not include the Gaussian factors of the wave 
functions Xf when constructing Eq. [^but only supple¬ 
ment the final result with the correct Gaussian factor. 
In addition, the wave function should be projected to 
the LLL for comparison with exact diagonalization re¬ 
sults. For the state described by Eq. [^with spin-up 
bosons and spin-down bosons, the two spin states 
have filling factors v"' = N^/[Nl{f^ + l)//t + and 
’'^ = Nt/[Nt{h-l)/h + Nl]. 

An alternative understanding of the trial wave func- 
tions can be obtained using the parton const ructiorP^Illl, 
In this framework, one physical boson is decomposed to 
multiple fictitious partons, which move in the same mag¬ 
netic field as the physical bosons and form IQH states. 
The state of the physical bosons is obtained by gluing the 
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FIG. 1. Energy spectra of bosons on the sphere. The sys¬ 
tem parameters are given as N^, A^, A^] on top of each 
panel and as [ff, fi] inside each panel. The dots and asterisks 
represent the energies of the trial wave functions Eq. In 
some cases, these two types of symbols have very close energy 
values and may not be easy to discern by inspection. 
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EIG. 2. Energy spectra of bosons on the torus. The aspect 
ratio L\/L 2 is one for all panels. The system parameters are 
given as [N^, A^, A^, A^] on top of each panel and as [/^, /j,] 
inside each panel. The number 2 in panel (b) indicates two 
quasi-degenerate states which can not be resolved by eye. 


partons together. When the partons form IQH states, 
we may obtain gapped topological phases of the physical 
bosons. This interpretation requires so we have 

to choose and to be different. For the trial wave 
function Eq. one physical boson with spin a consists 
of one type I parton with spin cr and one spinless type H 
parton. The spin-up (spin-down) type I parton form the 
^/t state and the spinless type H parton forms 

the xi state. 
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TABLE I. The overlaps between trial wave functions and exact eigenstates shown in Fig. L is the total angular momentum 
quantum number. The total number of linearly independent L multiplets is given in parentheses below each overlap. The 
entries with one (two) number(s) are the overlaps between the trial wave functions represented by the dots (asterisks) and the 
exact states. 


(a) [1,1] 


(b) [2,1] 


0.5 

0.4 


*■ 

* 

*r 


1 

9 

i 

♦ 

■ 

0.2 




0.1 











0.4 

* s 


* 

0.3 

■ T 1.0 

* 

1 

^ 0.2 

0.1 

• T 1.5 

♦ T2.0 

* 

♦ 


Nb = 

12 

8 


_u 

1/Nb 


FIG. 3. Energy gap values on the sphere and the torus. The 
horizontal axis is the inverse of the total number of bosons 
Nb = + N^. The values on the sphere have been rescaled 

by a factor to account for the deviation of the 

magnetic lengtlP^. (a) The /t = l)/ 4 . = 1 state; (b) the 
= 2, f_i = 1 state. The legend is shown in panel (b) where 
“S” means sphere, “T” means torus, and the number is the 
aspect ratio I/ 1 /I /2 of the torus. 


IV. NUMERICAL RESULTS 

We show in Fig. and Fig. the energy spectra of 
bosons on sphere and torus with four difference choices of 
ft and fi in Eq. 0 For bosons on the sphere, we are able 
to construct the trial wave functions explicitly. The over¬ 
laps between the trial states and the exact eigenstates are 
shown in Table |Tj The trial wave functions for the ground 
states (black squares in Fig. are excellent approxima¬ 
tions. The trial wave functions for the excited states 
can be obtained by promoting one composite fermion to 
the lowest empty composite fermion LL. For two compo¬ 
nent bosons in the same artihcial gauge potentiald^, one 
should consider all possible ways of exciting one compos¬ 
ite fermion to account for the low-lying excitations. The 
same construction in our system results in the blue aster¬ 
isks in Fig. (i.e. the states with one excited spin-up or 
spin-down composite fermion). The asterisks are accu¬ 
rate approximations of the low-lying excitations in Fig. 
n(c) and (d), but some of them have quite high energies 
^o do not represent the low-lying excitations) in Fig. 

(a) and (b). To further reveal the difference between our 
systems and those studied before, we construct the trial 
wave functions with only one spin-up composite fermion 
being excited and give their energies as red dots in Fig. 


which turn out to be a reasonably good description 
of the low-lying excitations in Fig. (a) and (b) but 
not in Fig. (c) and (d). In some angular momentum 
sectors, these two methods result in the same states so 
the dots and asterisks coincide with each other. The fact 
that the neutral excitations of the = 1,/^ = 1 and 
the = 2, fi = 1 states can be modeled very well using 
one type of composite fermions makes them very differ¬ 
ent from those studied in previous work^^. The energy 
gap values for several different systems corresponding to 
the = 1,/^ = 1 and the /^ = 2,/^ = 1 states are 
plotted in Fig. Although it is difficult to perform an 
accurate finite size scaling analysis, it is plausible that 
these two states are gapped in the thermodynamic limit. 
We also note that finite size effects are stronger on the 
sphere because the states appear to have larger gaps on 
the sphere. 

It is useful to consider the limit where and/or 
become very large. For the case with both —)■ 00 , the 

trial wave function Eq. actually describes composite 
fermions in zero effective magnetic field that form a gap¬ 
less Eermi liquid whose low-lying excitations contain ex¬ 
cited spin-up and spin-down composite fermion^^. If we 
choose one of to be small and the other to be large, it 
is very likely that the system would also be gapless due to 
the composite fermions with large effective filling factor. 
We expect that Fermi-liquid-like features will show up in 
a system with sufficiently large and/or /^. The precise 
values at which the system ceases to be gapped can not be 
determined for two reasons: (1) the gap value for a spe¬ 
cific set of and fi should be extracted from finite size 
scaling but there are not enough available data points in 
many cases; (2) the systems with large and/or can 
only be defined for a large number of bosons so can not 
be studied using exact diagonalization. We note that the 
energy spectrum of the = 2, fi = 2 state on the torus 
[Fig. I (d)] has a small energy gap. This is further sup¬ 
ported by effective field theory analysis presented below. 
It is possible that the /-^ = 2, /^ = 2 state may actually 
be gapless in the thermodynamic limit. The fact that 
its low-lying excitations are better modeled using com¬ 
posite fermions of both types is also consistent with the 
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behavior in the limit of infinite and f^. 

However, the analysis presented above does not seem 
to be applicable to the = 1, fi = 2 state. We can not 
claim that the system is gapless in the thermodynamic 
limit based on the numerical evidence. On the other 


hand, its low-lying excitations are different from those 
of the probably gapped states at = 1,/^ = 1 and 
= 2,/^ = 1, which suggests that a different physical 
interpretation for this state is required. The effective 
field theories presented below provide some hints about 
how to resolve this issue. 


V. EFFECTIVE FIELD THEORIES 


We construct effective field theories for our systems in two ways. In the first approach, we take the charges of 
all bosons to be unit and use the Chern-Simons transformation to implement flux attachmeniPsES. The partition 
function of the system is Z = J exp(j f dt d^r C) with Lagrangian density 




idt - — (p - 
‘ 2M ’ 


J d'^r'p„{t,r)V„r{r-r')pr{t,r') 


(9) 


CT.T=t4. 


where pa(t,r) = is the density of bosons with spin a. We define ^/>(r) = 

exp \—i f d^r' arg(r — r')p(r)] with arg(r — r') being the angle between the vector r — r' and the x axis. The 
Lagrangian density is transformed to 






(idt - ao) - — (p - A'" + a) 




^ ^ ^ r 

+ - - Y d^r'p^{t,r)Var{r - r')pr{t,r') 


( 10 ) 


CT,r=tt 


where pa{t,r) = tlj*{t,r)ipa{t,'r), the Chern-Simons vector potential a(t, r) = f d^r' p{t,v')\z x (r — r')]/|r — r'p, and 
the Chern-Simons scalar potential oq is a Lagrangian multiplier. 


In the first line of Eq. [TOj we introduce a weak electromagnetic potential ^ in addition to AZ to obtain 


^CF = Y 

cr=td 


{idt -ao + o) “ ^ (P “ + a - A^,,^ 


'da{t,v) 


( 11 ) 


In the second line of Eq. 10 the interaction term can be written as gp{t,v)5{v — r')p{t,r') with p = p^ + p^ = 
(27r)“^V X a so we have 

£cs = ^e\ti,^ax{t, Y)d^a^{t, ^ J [V x a(t, r)] 5(r - r') [V x a(t, r')] 

The current is related to Ag,;^ via the linear response equation and the magnetoplasmon 

excitations are related to the poles of the response function The partition function can be written as = 

f ZI[a] exp(zS'cs) f V"] exp(iS'cF) with ^cf = f d?r Ccf and S'cs = / d^r Ccs- The first step is to integrate 
out the composite fermion field '0* and 0 to obtain Z[Ag^^] = / D[a] exp(*S'eff). The second step is to integrate out 
the Chern-Simons gauge field in to obtain Z[Agxt] = exp(zS'ext)- The gap of the magnetoplasmon mode at zero 
momentum is (see Appendix B for more details) 


An — 


\J (/t + l)^!^! (•^I ~ + /t ~ A + l)wt‘^I ~ (/t + l)*^! ~ ifi ~ 1)^1 


( 12 ) 


where u)u is the effective cyclotron frequency of composite fermions. It should be emphasized that the gap values 
obtained in exact diagonalization generally are not at zero momentum. The finite momentum gap values are usually 
smaller than the zero momentum one s so the latter still provides useful information about whether the system is 


gapped. For the case with = fi^ Eq. 


12 


simplifies to Aq = {^J + 1 — /t)a; which suggests that the gap descreases 

as increases. This is consistent with^the numerical result that the gap at /t = !> /t = 1 is more robust than the 
gap at = 2,fi = 2. Another important observation is that the gap at = 1, fi = 2 vanishes, which suggests that 
the gap as seen in numerical calculations may be due to strong finite size effects. In general, Eq. [T^ suggests that the 
energy gaps for all the states at fi = + 1 are zero. 
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The ground state degeneracy of the many-body states on the torus can also be computed using the flux attachment 
approach together with the functional bosonization methocP^, but we defer this calculation to Appendix B because we 
find that it can be obtained more easily using the parton constructioiP^. The spin-up type I parton, spin-down type 
I parton, and spinless type II parton have charges ei^ = - f^), = -1/(2/; - 1), and en = 2/;/(2/; - I), 

respectively. The gauge potentials are the same for all the partons so we omit the superscript. We denote the 
parton fields as ipia and ipu to construct the Lagrangian density 


-^parton — ^ 


{idt - «o) - ^ (P - ei^A -h a) 






{idt + “o) “ ^ (P “ 




(13) 


where the gauge field a^{t,v) is used to glue the partons together to satisfy the constraint X)(T=td “ '/’n/’ii- 

After integrating out the parton fields, we find that the low-energy effective theory for the state is a Chern- 

Simons theory with action £eff = (/; — /t + . This predicts that the ground state degeneracy of 

the state on the torus is /; — /; + I, which agrees with the numerical results for the /; = !,/; = ! and 

/; = 2, /; = I states. The parton theory does not give information about the energetics so it can not tell us that 
the /; = 2,/; = 2 state has a small gap or may even be gapless. For the /; = 1,/; = 2 state, the coefficient in 
the Chern-Simons action vanishes, which suggests that it does not represent a gapped topological phase. One may 
compute to higher orders of the gauge field a^(t, r) and the resulting theory is likely to be a Maxwell theory describing 
a superfluid state. This is also the case for all other states at /; = /; -h 1. 



FIG. 4. A schematic phase diagram of our system. The pa¬ 
rameters are given as [/;,/;] on the figure. The system is 
gapless for sufficiently large /; and/or /; but can be gapped 
if these two parameters are small enough. The states with 
/i = /t + 1 rn^y represent a class of exotic superfiuid. The 
[/;, /;] = [2, 2] state is shown twice because it may be gapped 
or gapless. The question marks are used to indicate that the 
nature of some states is not fully settled. 


VI. CONCLUSION AND DISCUSSION 

In conclusion, we have studied possible topological 
phases of two-component bosons in artificial gauge po¬ 
tentials whose magnitudes for the two components are 
different. We propose trial wave functions using the com¬ 
posite fermion theory to describe systems with spin-up 
and spin-down composite fermions moving in effective 


magnetic fields pointing to opposite directions. An alter¬ 
native interpretation is given using the parton construc¬ 
tion. The trial wave functions are found to be excellent 
approximations of exact diagonalization results for sev¬ 
eral systems with contact interaction. An interesting ob¬ 
servation is that the low-lying excitations of some states 
can be modeled by exciting only the spin-up composite 
fermions. We obtain effective field theories using flux at¬ 
tachment and parton construction, which might seem dif¬ 
ferent but they both predict correct ground state degen¬ 
eracies as confirmed by numerical calculations. The flux 
attachment method also gives useful information about 
the energetics of the states and helps us to understand 
the numerical results. Our results are summarized in the 
schematic phase diagram Fig. The system is gapless 
for sufficiently large /; and/or /; but can be gapped if 
these two parameters are small enough. It is possible 
that the /; = 2, /; = 2 state is already gapless. The 
states with /; = /; -|- 1 may represent a class of exotic 
superfluid and deserve further investigations. 

Our theoretical work points out a direction for the 
experimental study of artificial gauge potentials. The 
synthetic magnetic field in a rotating Bose-Einstein con¬ 
densate depends on the angular velocity, so two internal 
states will experience different fields if they are rotated 
at different speeds. It is unclear whether this can be 
achieved without destabilizing the system. A more prob¬ 
able avenue is to combine two artificial gauge potentials 
with one being the same for all bosons while the other has 
opposite directions for the two components. The experi¬ 
mental evidences of the states studied here are similar to 
those of other quantum Hall states. The compressibility 
of a cold atom system can be measured to see whether it 
is gappecP^. The low-lying excitations can be probed us- 
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ing spectroscopic technique^. The gapless edge modes 
localized at the boundary of a gapped topological state 
can be visualized in experiments^. 

Finally, we would like to suggest a few directions that 
may be explored in future works. In this paper, we 
have only studied two-component bosons with contact 
interaction and the excitations of Eq. do not possess 
non-Abelian braid statistics. It is natural to ask what 
happens if the interaction potential between bosons is 
long-ranged and how to design the artificial gauge po¬ 
tentials to facilitate the emergence of non-Abelian topo¬ 
logical states. One may also study bosonic or fermionic 
systems with more internal degrees of freedom coupled 
to general Abelian or non-Abelian artificial gauge poten¬ 
tials. 
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Appendix A: Hamiltonian Matrix Elements 

Eor the spherical geometry, the matrix elements are given by 

where r = i?(sin 9 cos 0, sin 6 sin (j), cos 9) and = r/R. One can show that 

V'mi ~ ( 1) * ^rax,m2 4’mx+m'2 




{n; + + 1 ) 


1/2 


CG 


471 -h m 

where CG = {N^ /2, —mi; A^/2, —m 2 \N ^/2 + A^/2, —mi — m 2 ) are Clebsch-Gordon coefficients, so we have 


p(7TTa _ I 7Q 


,n: 


>N. 


Nl 






— 7-10*^7711,m2 *^m3,m4 

1 


d^li 


t. 


N'i+Nl 

mi+m2 


(Ol 


t. 


n;+n; 

m3+m4 


(fll 


Nl NZ 

— ^2 ^’^i+^2>m3-l-m4*^mi,m2 *^m3,m4 

Eor the torus geometry, the matrix elements are given by 


Tp(7TTa 

-^mim2m4m3 


= j d^Tid^V2 [V'^^(ri)]* [V';ijr2)]*(5(ri-r2)V'™^(r2)V'^3(ri) 


where r = {x,y) is the usual two-dimensional coordinates. We define the reciprocal lattice vectors as Gi = 27rea;/Li 
and G 2 = 27rey / L 2 ■ The interaction potential can be expressed in momentum space as 


’'<■■‘“‘■=' = 1711 If 


1112 


where q = 91 G 1 -I- 9262 . One can show that 


,iq-(ri-r 2 ) 




d^Y [iprmir)]* = exp <{ -W^'b \ exp <{ + m 3 ) ) sZ,Zm 2 +q 2 


so we have 


T^(7TT(7 

■^mim2m4m3 


where 5^ * and has been defined in the main text. 


= ^ 


>1 


exp < i2'Kqi 


Ns 


(mi - 92/2) (m2 -f 92/2) 


m 


Nl 


^mi-l-m 2 ,m 3 -|-m4 

















Appendix B: Electromagnetic Response Function 


In this section, we will use r to denote space-time coordinates (ro = t, ri = x, r 2 = y) and q to denote frequency- 
momentum variables (go = w, qi = Qxi 92 = %)■ The Chern-Simons gauge field a can be separated into an average 
(a) and a fluctuation a'^. The (a) part reduces the gauge field A'" for physical bosons to an effective gauge field A 
for composite fermions. We denote the efi'ective magnetic field strength for composite fermions as B , the effective 
cyclotron frequency of composite fermions as uja, and the filling factor of composite fermions as (they are the m 
and n used in the main text). We introduce b = a — Ag^^j and Si = J dS'r Cj with 


Cl = -b^ip*{r)i;a{r) - —h'^ ■ ^p*{r) (^p - A'") - (p + A'')^’*(r)V'^(r) 

The effective action SeS can be computed to the second order of Sj as 




SeS = Scs + (Si) + ^ (Sj) = Scs + l f ^ r2)6^(r2) 

where is the polarization tensor. In terms of the Green’s function of the composite fermion field 

Gi 2 = G'^{ri,r2) = -i (^^a{riWAr2)] 
the components of the polarization tensor are 
noo(^i)^2) = 'iGi2G2i 

K^irl,r2) = ^ [g?2(V2 - - (V^ + *A"),G?2G^i 

n-o(ri,r2) = ^ [(Vi - *A"),G?2 G^i - + *A"),G^i 

n-(rr,r2) = -^M(n-^2)-^ 

(Vi - zA"),G? 2 (V 2 - iA"),G^i + (V 2 + iA"),G? 2 (Vi + 

-(Vi - zA"),(V 2 + tA"),G^2G^i - GU^i + zA''),(V2 - lA^ )^G^, 


where Vi (i = 1, 2) takes derivative with respect to ri. 

It is convenient to transform the effective action to momentum space. The first part becomes 


f 

ao{q) 

/ 9 [ao(-<?),q-a(-g),qx a(-g)] G(g) 

q • a(g) 
q X a(g) 


where q is the unit momentum vector, q x a(g) denotes its 2 component (the only non-zero component), and the 
matrix G{q) is 

0 0 i0\q\ 

0 0 i9qo 

-i9\q\ -i9qo -0‘^\q\‘^V{q) 

with 9 = (27r)“^ and V{q) being the Fourier transform of the interaction. The second part becomes 


\ E b;{-qW^Aqmq) 

'J _•+• I 


<T=td 


and the components of the polarization tensor are 

Koiq) = \q\^FSiq) 

n«(9) = qoqzFo (g) - id^e.jqjF^ (g) 

^laiq) = qoq^Fo (g) -f idaitjqjF^ (g) 

n?,(g) = 5^3qlFo (g) + id^e^^q^F^ (g) -f ((5*j |gp - q^qJ)F^ (g) 
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where da denotes the direction of the effective magnetic field A . The three quantities FQ i 2 iQ) 

t\ 


Fo^q) = - 


\B |e 
27rM 


F^q) = 


2|B-| 


oo r-1 

■EE 

OO r -1 


s — t 


ql-(s-tYCols\ \2\B' 


L 


s-t I q 
2\B° 


s — t 


tl 


27rM2 


^ ^ q^-{s-t)^(Dls\ \2\B° 


s=f'^ t=0 


S — t— 1 


2\B 


q 

— (7 I 


2\B 


L 


2 

— a I 


F^^{q) = - 


\B |e 
27rM3 


(A. 

\2\B 

' oo r-1 

■EE 

s=f<’ t=0 


+ 2 L(rj+^ 

s — t 


2\B^\ 

t\ 


(1 - <5t,o) 


— (s — t)L 


s-t I q 
2\B° 


^s-t / q 


2\B 


s —i+1 


+ 2LtZl 


ql - {s-tyujl s! v 2 |b’ 
( q^ 


\2\B^ 


2\B 


L 


.-tf q^ 


\2\B° 


+ 2LlZl^^ 


(1 - ^qo) 

( q^ 


\2\B^ 


(1 - Ko) 


-2{s-t)Ll 


2\B 


where L^{x) is the associated Laguerre polynomial. 
We now have the effective action 


r 

ao{q) 

/ d 9 [ao(-9), q • a(-9), q x a(-g)] D{q) 

q • a(g) 

J 

q X a{q) 


“E/^^9[«o(-9),q-a(-g),q X a(-g)]r^(g) 


24ext.o(<7) 
q- Aext(9) 
q>^KKt{q) 


+ rxt,o(-9)>9'Aext(-9)>9 X Aext(-<?)] Ea{q) 


^Sxt.o(9) 

q-K^ti.q) 
9 X KAq) 


where 


Ea{q) 


Fo{q)\q? F^{q)qo\q\ idaF^{q)\q\ 

FE{q)qo\q\ ^o''(9)9o idaF[{q)qo 

-idaF[{q)\q\ -idaF[{q)qo F^{q)q'^ + F^{q)\q\‘^ 


and D{q) = C{q) + X]o'=-i' ^ Ea{q). The inverse of D{q) is 


D-\q) = 


X 


1 


{ql + m [F§{q)ql + Foiq)F,{q)\q\^ - F^q)] 

^ [Fo{q)q^O + E2{q)\q\^] [Fo{q)ql + ^ 2 ( 9 ) 19 ^] 


q2 + |q|2 ^ - J 52 + |^|2 

[-^ 0 ( 9)90 + E2{q)\q\'^] ^1^ + ^2(9)l9p] 

*-Fi(9)|9| *-^ 1 ( 9)90 


-*-F'i(9)|9| 
- 1 ^ 1 ( 9)90 
-^0(9) ( 9 o + | 9 P) - 


with Fo{q) = J2a=r,iE^{q), Fi{q) = ^ 2 ( 9 ) = Ea=td-^ 2 ^( 9 ) - 0‘^V{q). One can replace 

the Chern-Simons gauge field a in f D[a\ by its fluctuation a and the integration gives 


f d^q [A2xt,o(-9),9-X A2xt(-9)] df‘^^(9) 

1- 

q 

r 

L_ 


’ ^ 1 

<7T=tE 

_?x 


where the electromagnetic response function K'^^{q) = Ea{q)5aT —Ea{q)D ^{q)Er{q). The dispersion of the collective 
magnetoplasma modes is determined by the poles of K'^'^{q) at which 

F^{q)ql + F,{q)F2{q)\q\^-FKq)=^ 
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One can solve this equation and get the zero momentum gap 

\l (/t + ~ + 2(/tA + /t ~ A + ~ (/t +~ (A ~ 1)^1 

The ground state degeneracy of the many-body states on torus can also be obtained in this approach. In the small 
q limit, the action becomes 

^[Aext] = ^ J d^r 



with 


= d^rs^r 


dgdrrr 

EaA/‘" + l 


We introduce a gauge field b to express the particle current as /{2'k) and using the functional bosoniza- 

tion method to obtain the effective action 

*S'eff = -^ J d^r Q'^'^ex^ub'{{r)d^bl{r) 


with Q = G The ground state degeneracy of the state on torus is |det(5|. For the = !> A = 1 state, the 
Q-matrix is 



which predicts one ground state on torus. For the /^ = 2, A = 1 state, the Q-matrix is 

3/2 1 
1 0 


However, this matrix does not give a gauge-invariant Chern-Simons theory because = 3/2 is not an integer. To 

fix this problem, we introduce an auxilary gauge field 6 ^ and enlarge b to The effective action can be 

rewritten as 


with 




/2 1 0 
Q= \ 12 1 
\0 1 0 


which predicts two ground states on torus. 
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